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Abstract 

We discuss energy dependence of gap survival probability which follows 
from rational form of amplitude unitarization. In contrast to eikonal form of 
unitarization which leads to decreasing energy dependence of gap survival 
probability, we predict a non-monotonous form for this dependence. 
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Studies of the processes with large rapidity gaps are very important as a tool 
for the search of new physics. For the first time such processes have been dis- 
cussed in m 121 13. Predictive power of QCD calculations for cross-sections 
of such processes is affected by the uncertainties related to the soft interactions 
(rescattering) in initial and final states. Dynamics of such interactions is accounted 
then by the introduction of a factor which is known as a gap survival probability 
0, i.e. a probability to keep away inelastic interactions which can result in filling 
up by hadrons the large rapidity gaps. Energy dependence and magnitude of gap 
survival probability is an important issue e.g. in the studies of Higgs production in 
double diffractive exclusive and inclusive processes at Tevatron and the LHC (cf. 

The extensive studies of magnitude and energy dependence of gap survival 
probability have been performed, results of these studies can be found, e.g. in 
papers fSHIIllHl. 

The gap survival probability {\S\'^) is determined by the relation Q: 

.^.,, Io^DH{b)\S{s,b)\'d% 

where Dnib) is the probability to observe a specific hard interaction in collision 
of the hadrons hi and /i2 and P{s, b) = \S{s, where S is elastic scattering 
^-matrix, i.e. P is a probability of the absence of the inelastic interactions. In the 
eikonal formalism which is usually used for estimation of {\S\'^) the probability 
P{s,b) = exp(— 6)). All estimations of the gap survival probability per- 
formed on the basis of the eikonal amplitude unitarization lead to decreasing en- 
ergy dependence of this quantity. Therefore rather small values of cross-sections 
for diffractive Higgs productions are expected at the LHC energies ll9l fT0ll . 

However, there is an alternative approach to unitarization which utilizes a ra- 
tional representation and leads, as it will be shown below, to a non-monotonous 
energy dependence of gap survival probability. Arguments based on analytical 
properties of the scattering amplitude ifTTIl provide support for the rational form 
of unitarization. In potential scattering rational form of unitarization corresponds 
to an approximate wave function which changes both the phase and amplitude of 
the wave. The rational form of unitarization in quantum field theory is based on 
the relativistic generalization lfT2ll of the Heitler equation llT3l . In the [/-matrix 
approach based on rational form of unitarization, the elastic scattering amplitude 
in the impact parameter representation has the form: 

= (2) 

1 — zU{s, b) 
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where U{s, h) is the generalized reaction matrix, which is considered to be an 
input dynamical quantity similar to the eikonal function. Unitarity equation for 
the elastic amplitude f{s,b) rewritten at high energies has the form 

Imf{s,b) = \f{s,b)\' + v{s,b) (3) 

where the inelastic overlap function 

1 d^inel 



r]{s,b) 



is the sum of all inelastic channel contributions. Inelastic overlap function is re- 
lated to U{s, b) according to Eqs. ^ and Q as follows 

vis b) = ^"^^^^'^^ (4) 



The probability 

l + iU{s,b) 



Pis,b)^\Sis,b)\' 



tU{s,b) 



2 

(5) 



Unitarity of the scattering matrix implies, in principle, an existence at high enough 
energies s > sq, where sq is the threshold, of the new scattering mode - anti- 
shadow one. It has been revealed in [ 14| and effects related to antishadowing at 
the LHC energies have been discussed in [15|. The most important feature of 
this mode is the self-damping of the contribution from the inelastic channels. The 
rational form of unitarization provides smooth transition beyond the black disk 
limit, where antishadow scattering mode is realized, i.e. at high energies and at 
small impact parameters elastic scattering channel can play dominating role. 

There is an experimental indication that this mode can indeed occur at very 
high energies. The analysis of the experimental data on high-energy diffractive 
scattering shows that the effective interaction area expands with energy and the 
interaction intensity - opacity - increases with energy at fixed impact parameter. 
At the Tevatron highest energy ^/s = 1800 GeV elastic scattering amplitude is 
very close to black disk limit at 6 = L16J . i.e. 

Im/(s, 6 = 0) = 0.492 ± 0.008. 

The central inelastic collisions of hadrons are far from amalgamation of the two 
hadrons in one region of space as it was shown in IfTTIl and the persistence of 
longitudinal momentum takes place at very high energies. 
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The rational form of amplitude unitarization can be put into agreement with 
the experimental data using various model parameterizations for the [/-matrix. 
They all lead to the same qualitative predictions, which reflect general properties 
of this unitarization scheme. Originally the Regge-pole model was used to get 
an explicit form of fZ-matrix [ 18 J and a good description of the experimental data 
has been obtained in this model ifT^ as well as in its Dipole-Pomeron modification 

mi- 

We use the model for the tZ-matrix based on the ideas of chiral quark ap- 
proaches f2l\. It is in a good agreement with the data 1.15. .22. .23 1 and is also 
applicable to the large angle scattering. We would like to stress here that the qual- 
itative conclusions of the present paper do not depend on the particular tZ-matrix 
parameterization. 

In this model the picture of a hadron consisting of constituent quarks embed- 
ded into quark condensate is used. This picture implies that overlapping and inter- 
action of peripheral clouds occur at the first stage of hadron interaction. Nonlinear 
field couplings could transform then the kinetic energy to internal energy and as 
a result massive virtual quarks appear in the overlapping region. These quarks 
generate an effective field. Valence constituent quarks located in the central part 
of hadrons are supposed to scatter simultaneously in a quasi-independent way in 
this field. Massive virtual quarks play a role of scatterers for the valence quarks 
in elastic scattering and their hadronization leads to soft production process of 
secondary particles in the central region [24|. The number of such scatterers was 
estimated 

N(s,b) 0. (i^iMVfBcW, (6) 

under assumption that part of hadron energy carried by the outer condensate 
clouds released in the overlap region to generate massive quarks, where mg is 
a constituent quark mass and (/cg) - an average fraction of hadron energy carried 
by the constituent valence quarks. Function Dc{h) is a convolution of the two 
condensate distributions D^^{h) and D^'^{h) inside the hadron hi and /i2. 

We will consider for simplicity the case of a pure imaginary amplitude, i.e. 
U iu. The function u{s,b) is represented in the model as a product of the 
averaged quark amplitudes (/q), 

N 

u{s,b) = l[{fQXs,b)) (7) 

i=l 

in accordance with assumed quasi-independent nature of the valence quark scat- 
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tering, is the total number of valence quarks in the colliding hadrons. The 
essential point here is the rise with energy of the number of the scatterers like y/s. 
The 6-dependence of the function (/q) has a simple form (/q(6)) oc exp(— rngfe/^). 
The generalized reaction matrix gets the following form 



u{s, b) = g i 1 + a 



N 



exp( 



Mb, 
'T' 



(8) 



where M = J2q=i "^Q- Here niq is the mass of constituent quark, which is taken 
to be 0.35 GeV\ 

This model provides linear dependence on ^/s for the total cross-sections, i.e. 
<^tot = a + Cy/s in the limited energy range ^/s < 0.5 TeV. Such behaviour and 
model predictions for higher energies are as it was already mentioned in agree- 
ment (Fig. 1) with the experimental data on total, elastic and diffractive scattering 
cross-sections lE^E^ . 





10' W 10= 10* s"^GeV 



10' S"^GeV 



Figure 1 : Total and ratio of elastic to total cross-sections of pp and pp-interactions 

This unitarization approach leads to the following asymptotical dependencies 
a tot oc In^ s and ainei oc In s, which are the same for the various models and reflect 
essential properties of this unitarization scheme. 

Thus, now the probability P(s, b) = \S{s, 6)p can be calculated in a straight- 
forward way, i.e we use for the function u{s, b) formula (HJ) with parameter fixed 
from the total cross section fit and relation of S{s, b) and U -matrix (jSj). The impact 
parameter dependence of P(s, b) for the different energies is presented on Fig. 2. 

'Other parameters have the following values: g = 0.24, ^ = 2.5, a = 0.56 • 10~^ which have 
been obtained from the fit to the total hp cross sections II22I . 
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Figure 2: Impact parameter dependence of probability P{s,b) = |5(s,6)p at three 
values of energy ^/s = 500 GeV (shadow scattering mode), y/s = 1800 GeV (black disk 
limit at b = 0) and ^/s = 14 TeV (shadow and antishadow scattering modes). 

To calculate the gap survival probability (IS'p) we need to know probability of 
hard interactions Dnib). To be specific we consider the hard central production 
processes 

p + p ^ p + gap + {Higgs or jj) + gap + p (9) 

The interest in such processes is related to the clear experimental signature and 
significant signal-to-background ratio. 

We can write down the probability of hard interactions in the model as a con- 
volution 

Dnib) = aH j D^^{h^)wH{h + bi - b2)Z)^(b2)rfbirfb2, (10) 

where ^^^^(b + bi — b2) is the probability of hard condensate (parton) interactions. 
It is natural to assume that the probability wh has much steeper impact parameter 
dependence than the functions D^^ (b) have and, therefore, the impact parameter 
dependence of wh determines behaviour of Dnib). Thus, we assume a simple 
exponential dependence for the function Dni^b), i.e. 

DH{b)-aHexp{-MHb), (11) 

where mass Mh is determined by the hard scale of the process. We perform 
numerical calculations of the gap survival probability (15^) using Eqs. ([T]), ^ 
and (Hill- 
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We take the hard scale to be determined by the mass of J/\l/-meson, i.e. 
Mh = Mj/y^,. This choice is in accord with the fact that production of J/\l/- 
meson can be treated as a hard process and therefore its mass set a hard scale ll25l . 
Lower values of Mh lie in the soft region ^. It should be noted that numerical 
results are rather stable and depend weakly on the scale Mh, when it is in the 
hard region, i.e. Mh > Mj/,$ e.g. for illustration we used the value Mh = 8 
GeV and it leads to slightly lower values of the gap survival probability at low 
energies. Results of calculations are presented in Fig. 3. One can notice that the 
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Figure 3: Energy dependence of gap survival probability {\S\' 



gap survival probability reaches its minimal values at the Tevatron highest energy. 
It is not surprising since the scattering at this energy is very close to the black disk 
limit. 

Asymptotical behaviour of the gap survival probability has the form: 



{\S\') 



^Mh 



2mc 



Ins. 



(12) 



The two unitarization schemes ([/-matrix and eikonal) lead to different predic- 
tions for the gap survival probability in the limit s — > oo; eikonal unitarization 
predicts (ISp) = at s — > oo, while [/-matrix formalism gives (l^p) = 1. 
Latter is a result of transition to the antishadow scattering mode in the [/-matrix 
unitarization [14|, when the amplitude becomes \f{s,b)\ > 1/2 (in the case of 

^In Eq. (CB <7H can be interpreted as a probability of the hot spot formation under condensate 
interaction and Rh — 1/Mh is the radius of this hot spot. 
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imaginary eikonal the scattering amplitude never exceeds the black disk limit 
\f{s,b)\ < 1/2). It should be noted that the Froissart-Martin bound implies 
unitarity (not black disk) limit for the partial amplitudes. When the amplitude 
exceeds the black disk limit (in central collisions at high energies) then the scat- 
tering at such impact parameters turns out to be of an antishadow nature. In this 
antishadow scattering mode the elastic amplitude increases with decrease of the 
inelastic channels contribution at small impact parameters and most of inelastic 
interactions occur in the peripheral region; the inelastic overlap function has a pe- 
ripheral impact parameter profile, which is a main reason of the large gap survival 
probability. 

Numerical predictions for the gap survival probability obtained here depend 
on the particular parameterization for U -matrix, but the qualitative picture of the 
energy behaviour of (l^p) reflects transition to the new scattering mode at the 
LHC energies and is valid for various [/-matrix parameterizations which provide 
rising total cross-sections. One should note that the numerical values of (15^) at 
the Tevatron energies are in a qualitative agreement with a number of quantitative 
calculations performed in the eikonal approaches (cf. iHI). 

In the sense of gap survival probability situation should be more favorable at 
the LHC energies since the obtained numerical values of {\S\'^) at these energies 
are close to unity and this should lead to much higher cross-sections (by factor 
of 40 compare to the calculations based on the gap survival probability estima- 
tions in the framework of the eikonal model lE^ ) e.g. for Higgs production in 
double diffractive processes compared to the values obtained with eikonal based 
estimations of the gap survival probability. 

Thus, the antishadowing which appearance is expected at the LHC energies, 
should be correlated with enhancement of Higgs production cross-section in dou- 
ble diffractive scattering and this would significantly help in detecting of Higgs 
boson. 
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